Análise Matemática IV 

Problemas para as Aulas Práticas 



6 de Junho de 2005 



Semana 9 



1. Seja 
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Determine e íA e resolva o problema de valor inicial x' = Ax, x(l) = 

[ 

Resolução: 

Os valores próprios de A são as soluções da equação 

det(A - XI) = & ((5 - A) (3 - A) + l) (4 - A) = O A = 4 

Para calcular o(s) vector(es) próprio(s) associados, há que determinar v = (a,b,c) G 
tal que 



(A - 4/)v = <* 
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^ a - b = 



Tem-se então que 

v = (a, b, c) = (a, a, c) = (a, a, 0) + (0, 0, c) = a(l, 1, 0) + c(0, 0, 1) 

donde podemos escolher 

vi = (1,1,0) , v 2 = (0,0,1) 

Concluimos que a matriz A não é diagonalizável, mas é semelhante a uma matriz de 
Jordan com dois blocos, isto é 

A = SJS" 1 

em que, por exemplo 
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S = [ V! V g V 2 ] 

em que v g é um vector próprio generalizado. É então solução da equação 



(A - 4I)v = vi <^ 
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<í=> a - b = 1 



Tem-se então que 



(a, b, c) = (1 + b, b, c) = (1, 0, 0) + (b, b, 0) + (0, 0, c) = (1, 0, 0) + 6(1, 1, 0) + c(0, 0, 1) 



donde podemos escolher 
Tem-se então que 



S = 



1 1 

1 o o 

1 



«=(1,0,0) 

e S- x = 



1 

1 -1 o 
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Finalmente 

e tA = e íSJS-i = Se íJ S -l = 

ou seja 
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e tA = é l 



í+1 -t 
t 1-í 

1 



Assim sendo, a solução do PVI, x' = Ax, x(l) = 



1 






é dada por 



x(t) = e A ( í " 1 )x(l) = e 4 ^" 1 ) 



2. Seja 
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Resolva o problema de valor inicial x' = Ax, x(0) = 



Resolução: 



Visto A ser uma matriz triangular, os seus valores próprios são os elementos da sua 
diagonal principal, ou seja 



det(A - XI) = & A = 2 ou A = ou A = 1 



Análise Matemática IV 



Para calcular o vector próprio associado ao valor próprio A = 2, há que determinar 
v = (a, 6, c) GR 3 \ {(0, 0, 0)} tal que 



(A - 2/)v = O 
Tem-se então que 
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O b=c= 



donde podemos escolher 



v= (a,b,c) = (a, 0,0) = a(l, 0, 0) 



vi = (1,0,0) 



Para calcular o vector próprio associado ao valor próprio A = 0, há que determinar 
v = (a, b, c) GR 3 \ {(0, 0, 0)} tal que 



(A - 0/)v = e> 
Tem-se então que 



2-2 2 

1 
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-£4> a = b e c = 



donde podemos escolher 



v = (a, b, c) = (a, a, 0) = a(l, 1, 0) 



V2 = (1,1,0) 



Para calcular o vector próprio associado ao valor próprio A = 1, há que determinar 
v = (a, b, c) GR 3 \ {(0, 0, 0)} tal que 



(A-/)v = & 

Tem-se então que 
donde podemos escolher 



1 -2 2 
-1 1 
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<í=>6 = cea = 



v = (a,b,c) = (0,6,6) = 6(0,1,1) 



V3= (0,1,1) 

Concluimos que a matriz A é diagonalizável, ou seja 

A = SDS" 1 



em que, 



D 



2 

1 



S = [ vi v 2 v 3 ] = 



1 1 
1 1 
1 
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Finalmente 



e íA = gíSDS- 1 



= Se^S" 1 = 



ou seja 
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Assim sendo, a solução do PVI, x' = Ax, x(0) = 



x(t) = e A(í - o) x(0) = 



,2* „2t 



e"- —e lt + I e 
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é dada por 



e*-l 
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3. Seja 



2 1 

2 

-10 3 



Resolva o problema de valor inicial x' = Ax, x(0) = 



Resolução: 

Os valores próprios de A são as soluções da equação 

det(A - XI) = <* (2 - A) 2 (3 - A) = O A = 2ouA = 3 

Para calcular o vector próprio associado ao valor próprio A = 2, há que determinar 
v = (a, b, c) GR 3 \ {(0, 0, 0)} tal que 



(A - 2/)v = O 
Tem-se então que 
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-1 1 
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<í=> a = ceò = 



donde podemos escolher 



v = (a, b, c) = (a, 0, a) = a(l, 0, 1) 



vi = (1,0,1) 



Para calcular o vector próprio associado ao valor próprio A = 3, há que determinar 
v = (a, b, c) el 3 \ {(0, 0, 0)} tal que 



(A - 3/)v = O 
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Tem-se então que 



donde podemos escolher 



v = (a,b,c) = (0,0, c) = c(0,0,l) 



v 2 = (0,0,1) 



Concluímos que a matriz A não é diagonalizável mas é semelhante a uma matriz de 
Jordan com 2 blocos (o número de vectores próprios linearmente independentes), ou seja 



A = SJS" 1 



em que, por exemplo, 
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S = [ vi vf v 2 ] 

onde V! e v 2 são os vectores próprios atrás calculados e vf é um vector próprio ge- 
neralizado correspondente ao valor próprio A = 2. Para o calcular há que determinar 
v = (a, b, c) GR 3 \ {(0, 0, 0)} tal que 



(A - 2I)v = vi & 

Tem-se então que 
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<í=> c = a + 1 e 6=1 



v = (a, b, c) = (a, 1, a + 1) = a(l, 0, 1) + (0, 1, 1) 
donde podemos escolher 



Tem-se então que 



S = 



1 o 

1 o 

1 1 1 



vf= (0,1,1) 



s-^ 



1 o o 

1 o 

-1 -1 1 



Finalmente 



D tA = e tsjs-i = Se a s -i = 



e 



2í 



te 2t 

,2t 




e 2t 

e 2í _ e 3í (j. + ^ e 2í _ e 3í e 3í 



Assim sendo, a solução do PVI, x' = Ax, x(0) = 



x(í) = e A ( í -°)x(0) = 



-,2í 



te 



2í 




e 2t 

e 2í _ e 3í (j. + ^ e 2í _ e 3í e 3í 



é dada por 



e 2í (í+l) 



-,2í 



2í ^3í 



(í + 2)e 2í -e 
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4. Determine a solução geral do seguinte sistema de equações diferenciais: 



x' = Ux - 10y+l 
y' = 10x-2y + 2 



Sugestão: Determine primeiro uma solução particular constante. 
Resolução: 

Escrevendo o sistema na forma matricial 



X 




14 


-10 




X 
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10 
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2 



AX + C(í) 



(1) 



Dado que a equação é linear, uma sua solução será da forma 

X = Xh + X p 

em que Xh é a solução geral do problema homogéneo associado X' = AX, e X p é uma 
solução particular da equação X' = AX + C. Seguindo a sugestão, vamos procurar 
uma solução particular constante, ou seja, vamos procurar X p = (ci,c 2 ) que verifique a 
equação (1), isto é 



Cl 

c 2 



14 -10 
10 -2 





Cl 

c 2 


+ 


1 

2 



, 14ci- 10c 2 + l = , 11, 

* < 10ci-2c 2 + 2 = * (^^ = (-4' "4' 



pelo que X p = (— |, — |). De seguida calcularemos Xh, que como já sabemos envolve o 
cálculo da matriz e íA . Os valores próprios de A são as soluções da equação 

det[A - XI] = <* A 2 - 12A + 72 = O A = 6 ± 6« 

Para calcular o vector próprio associado ao valor próprio A = 6 + 6i, há que determinar 
v = (a, 6) 7^ (0, 0) tal que 



(A- (6 + 6í)/)v = O 
Tem-se então que 



6?' 



10 



(a, 6) = (a, 



-10 
-8-6i 

-3i , 



a 

b 



O 6 = 



3?; 



a(l, 



3r 



donde podemos escolher 



vi = (1, 



3r 



e como consequência, o vector próprio associado ao valor próprio A = 6 — 6i será 

3i, 



V2=(l, 



3i, 



(1- 



5 ' K ' 5 
Tem-se então, que A é uma matriz diagonalizável, isto é 

A = SDS" 1 
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em que 



pelo que 



Finalmente 



e tA = Se^S" 1 



D 



6 + 6« 
6-6i 



s- x = 



e S = 
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4-3i 4+3i 
5 5 



4 + 3« -5 
-4 + 3« 5 
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5 5 
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Q e (6-6i)í 



4 + 3í -5 

-4 + 3« 5 



Cif 
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6íf _ e -6ií) + 3i ( e 6ii + e -6ÍÍN 



25(-e 6íí + e- 6lí ) 



6it „—6it\ 



25 (e 



4sen(6í) + 3 cos(6í) — 5sen(6£) 



òU(-e 6it + e- 6it ) 



3i(e 



6ff i ~ — 6ít\ 



5sen(6í) — 4sen(6í) + 3 cos(6í) 



Temos então que 



M 



íAi 



X h = é A K = 



4sen(6í) + 3cos(6í) 



-5 sen(6í) 



5sen(6í) — 4sen(6í) + 3 cos(6í) 






e por fim a solução pedida é 



„6í 



X — Xh + Xp 



4sen(6í) + 3 cos(6í) — 5sen(6í) 



5sen(6í) 



-4sen(6í) + 3cos(6í) 






para k\, ki G HL 
5. Considere a seguinte matriz: 



1 1 
1 
1 



(a) Calcule e At . 
Resolução: 



,Aí 
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té 


0~ 
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(b) Determine a solução do problema de valor inicial 

y = Ay + h(t) 

. y(i) = (i,i,i) T 

onde h(t) = (0, 2e*, é) T . 
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Resolução: 

A solução geral da equação homogénea é dada por y — Ay = 



y(*) = e 



Aí 



(li 

«3 



Para determinar uma solução particular da equação basta determinar soluções par- 
ticulares das equações 

yi-yi-y2 = 

ÍJ3-y3 = e*. 

Segue-se que 

y 2 (t) = Até e y 3 (t) = Bté. 

Para determinar A e B tem-se 

m - y2 = 2e* => Aé + Até - Até = 2é => A = 2 



2/3 - 2/3 = é => 5e* + Síe* - 5íe* = e* => 5 = 1 
Substituindo y 2 = 2íe* na primeira equação, obtém-se: 

yi-yi = 2té 

Procedendo como anteriormente, resulta que y\(t) = t 2 e l é uma solução da equação 
acima. Logo uma solução particular é: 

~t 2 é~ 
y P (t)= 2te t 
té 

Aplicando a condição inicial, obtém-se 

y(í) = e A(i - 1J 

Nota: O problema também pode ser resolvido usando a fórmula da variação das 
constantes. 

6. Considere o seguinte problema de valor inicial 

y' = Ay + b(í) 



1 




"tV" 


l-2e 


+ 


2té 


1 -e 




_té_ 



y(0)= [0 ] : 



onde 



-2000 
3-2 
12 
0-2 1 



e b(í) = 
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(i) Determine a solução geral da equação homogénea. 

(ii) Sendo y(í) = [yi(t) y 2 (£) ys(t) Vii^Y a solução do problema não homogéneo, deter- 
mine y 2 (3). 

Resolução: 

(i) Visto a matriz A ser da forma 



Ai 
A 2 



sendo A x = 



-2 
3 -2 



e Ao = 



1 2 
-2 1 



temos que 



,Aí 



e Alí 
e A2Í 



o que facilita bastante os cálculos. Comecemos por calcular e Alt . O valor próprio de Ai 
é —2 associado aovector próprio vi = (0, 1) e ao vector próprio generalizado vf = (|, 0). 
Sendo assim 



Ai = SJS" 1 = 
e consequentemente 

e Alí = Se^S" 1 = 



I 
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1 o 



-2í 



te 







-2 1 
-2 



-2í 
-2í 
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" 1 " 
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= e 


u 3í 1 



Vamos agora calcular e A2Í . Os valores próprios de A 2 são 1 + 2i e 1 — 2i associados 
(respectivamente) aos vectores próprios v x = (l,i) e v 2 = (1, — i). Sendo assim 



A 2 = SDS" 1 = 

e consequentemente 



,A 3 í 



= Se^S" 1 = 



2?; 



1 1 

i —i 



1 1 
i —i 



:>it 



l + 2i 




Finalmente 




e~ 2it 



-2í 




l-2i 



i 1 
i -1 



1 
2i 



i 1 
i -1 



= e 



cos(2í) sen(2í) 
- sen(2í) cos(2í) 



,Aí 



,Aií 



e 



A 2 í 



3te" 2í e~ 2t 







e^os^í) e*sen(2í) 

-e í sen(2í) e*cos(2í) 



e a solução geral do sistema homogéneo é dada por 

-2* o 



y(í) = e AÍ C = 



■,Aií 







e 



A 2 í 



e - u 

3íe" 2í e" 2í 

e^os^í) e í sen(2í) 

-e í sen(2í) e*cos(2t) 



Cl 

c 2 

C 4 
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(ii) Pela fórmula da variação das constantes, a solução do problema de valor inicial dado 
será 



y(*) = e 



Aí 



o 
o 
o 

o 

-2* 



,Aí 



-Asi 



f "*" / e " s b(s) ds 
'o 







3te" 2í e" 2í 



r 



2e 2s 


7o 







rf.s 







e*cos(2í) e*sen(2í) 

-e*sen(2í) e*cos(2í) 



1 - e" 2í 







pelo que 

y 2 (3) = 1 - e" 6 

7. (i) Determine a solução do sistema linear 

x' = x — y 
y' = 2x - y 

que satisfaz a condição inicial x(0) = y(0) + 1 = 1. 
(ii) Considerando agora o sistema 

x 1 = x — y 

y' = 2x-y 

z' = y — (sent)z 

utilize a alínea anterior para determinar a solução que verifica a condição inicial 
x(0) = y(0) + l = z(0) = 1. 

Resolução: 

(i) Escrevendo o sistema na forma matricial 



x 



1 -1 

2 -1 



x = Ax 



Os valores próprios de A são i e —i associados (respectivamente) aos vectores próprios 
vi = (1, 1 — i) e V2 = (1, 1 + i). Sendo assim 



A = SDS 



2?; 



í í 

\-i 1 + i 



i 
-i 



1 + i -1 
-1 + i 1 
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e consequentemente 

e At = Se^S" 1 



1 
2i 



1 1 

1 - i 1 + i 



e u 
e~ H 



1 + i -1 
-1 + i 1 



sent + cosi — sení 

2sení — sen £ + cosi 



x(t) ' 


- P Kt 


" 1 " 




sení + cosi 


y(t) . 


— e 







2 sení 



pelo que 



(ii) Dado que nas duas primeiras equações não há dependência em z podemos concluir 
de imediato por (i) que 

x(t) = sení + cosi e y(t) = 2 sen. t 

Falta então determinar z, ou seja resolver o PVI 

z' = 2sent — (sent)z , z(0) = 1 

Trata-se de uma equação linear, de factor integrante 

/x(í) = e f sentdt = e~ cost 

Então, a equação é equivalente a 

': ) - 2 sen ff - cos< ^ e" cosí 2 = 2e- cosí + c ^ z(t) = 2 + ce 
Dado que z(0) = 1, conclui-se 



d / 

I — cost „ \ o ^^-^-/-^— cosi 



cosi 



z(i) = 2 - e 



cosi— 1 



